Anomalous Andreev bound state in non-centrosymmetric superconductors 
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We study edge states of non-centrosymmetric superconductors where spin-singlet d-wave pairing 
mixes with spin-triplet p (or /)-wave one by spin-orbit coupling. For dxy-wave pairing, the obtained 
Andreev bound state has an anomalous dispersion as compared to conventional helical edge modes. 
A unique topologically protected time-reversal invariant Majorana bound state appears at the edge. 
The charge conductance in the non-centrosymmetric superconductor junctions reflects the anoma- 
lous structures of the dispersions, particularly the time-reversal invariant Majorana bound state is 
manifested as a zero bias conductance peak. 

PACS numbers: 74.45. -fc, 74.50. +r, 74.20. Rp 
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Recently, physics of non-centrosymmetric (NCS) su- 
perconductors is one of the important issues in condensed 
matter physics [H-Q. One of the remarkable features in 
NCS superconductors is that due to the broken inver- 
sion symmetry, superconducting pair potential becomes 
a mixture of spin-singlet even-parity and spin-triplet odd- 
parity [1]. Due to the mixture of spin-singlet and spin- 
triplet pairings, several novel properties such as the large 
upper critical field are expected 0, 0] . 

Up to now, there have been several studies about su- 
perconducting profiles of NCS superconductors [sI-Ig- 13|. 
In these works, pairing symmetry of NCS superconduc- 
tors has been mainly assumed to be s -I- p-wave. How- 
ever, in a strongly correlated system, this assumption is 
not valid anymore. Microscopic calculations have shown 
that dj.2_y2-waYe spin-singlet pairing mixes with /-wave 
pairing based on the Hubbard model near half filling Id'] . 
Also, a possible pairing symmetry of superconductivity 
generated at heterointerface LaAlO^SrTiOs Q has been 
studied based on a similar model 15]. It has been found 
that the gap function consists of spin-singlet c?a;j,-wave 
component and spin-triplet p-wave one [15|. Therefore, 
now, it is a challenging issue to reveal novel properties 
specific to dxy -I- p or d^^-yi. -\- /-wave pairing. 

The generation of Andreev bound state(ABS) at the 
surface or interface is a significantly important feature 
specific to unconventional pairing since ABS directly 
manifests itself in the tunneling spectroscopy. Actually, 
for dxy-wave pairing, zero energy dispersionless ABS ap- 
pears [16]. The presence of ABS has been verified by 
tunneling experiments of high-Tc cuprate 17| as a zero 
bias conductance peak (ZBCP). For NCS superconduc- 
tors, when p-wave pair potential is larger than s-wave 
one, it has been shown that ABS is generated at its edge 
as helical edge modes similar to those in quantum spin 
Hall system [8 



12|. Several new features of spin trans- 



NCS superconductors. Since tunneling spectroscopy via 
ABS 16| is a powerful method to identify pairing sym- 
metry and mechanism of unconventional superconductors 
17| . it is quite important and interesting to clarify ABS 
and resulting tunneling conductance for dxy+p-wave and 
dx2^y2 -\- /-wave pairings. 

In this Letter, we investigate ABS and tunneling con- 
ductance ac in normal metal / NCS superconductor 
junctions. For both d^y -f- p-wave and dx2-y2 + /-wave 
cases, new types of ABS are obtained. In particular, 
for dxy + p-wave case, due to the Fermi surface split- 
ting by spin-orbit coupling, a single branch of topologi- 
cally stable Majorana bound state appears. Recently, to 
search for Majorana fermions is one of the hottest issues 
in condensed matter physics 3, 1^ • In stark contrast to 
the other Majorana fermions, the present one preserves 
time-reversal symmetry. From this difference, the "time- 
reversal invariant (TRI) Majorana bound state" has a 
peculiar flat dispersion. It shows a unique ZBCP in ac 
depending on the spin-orbit coupling. Therefore, the ex- 
perimental identification is feasible. 

We start with the Hamiltonian of NCS superconductor 
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port stemming from these helical edge modes have been 
also predicted However, there has been no the- 



ory about ABS in dxy + p- or dx2-y2 -\- /-wave pairing in 



with H{k) = Ck + V{k) ■ &, V{k) = \{xky - ykx), 
= n'^k'^/{2m) - /i. Here, /i, m, (T and A denote 
chemical potential, effective mass, Pauli matrices and 
coupling constant of Rashba spin-orbit interaction, re- 
spectively 0. The pair potential A(fc) is given by 
A(fc) — [d{k) ■ &]iay -\- iip{k)ay. Due to the spin-orbit 
coupling, the spin-triplet component d{k) is aligned with 
the polarization vector of the Rashba spin orbit coupling, 
(i(fc)||V(fe)Q. Then, the triplet component is d{k) = 
Atf{k){xky — ykx)/k with k — \/k^ while singlet com- 
ponent reads ipik) = Asf{k) with Aj > and As > 0. 
/(fe) is given by /(fc) = 2kxky/k'^ for dxv + p-wave and 
/(fc) = (kl - kl)/k'^ for dx2_y2 + /-wave.Q The super- 
conducting gaps are Ai = |Ai(fc)| and A2 = |A2(fc)| for 
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the two spin-split band with Ai(A;) = (At + As)/(fc) and 
A2(fc) = (At — As)/(fe), respectively, in homogeneous 
state [i]. 

Let us consider a wave function including ABS local- 
ized at the surface. Consider a two-dimensional semi- 
infinite superconductor on a; > where the surface is 
located at a: = 0. The corresponding wave function is 
given by ^Vi\ 

*s(a;) = [citpt exp{iq^^x) + c^tl)^ exp(-i(7f^a;) 
+ct^t exp(i(7^^x) -I- c^V2' exp(-ig2;,a;)] exp{ikyy), (2) 

fcl(2) 



we look at the At > As case. For A^ = with /? = 0, 
Eb — icky with some constant c for small ky (curve a in 
Fig. 1 (A)) as shown in the case of s-l-p-wave pairing 
since 77 = is satisfied. This type of ABS is called helical 
edge mode [ll, 12, 21 1. However, this condition is satis- 
fied only for Ag = and (3 — 0. In fact, Sb near ky = 
becomes absent in general as shown in curves a in Figs. 
1(B), (D) and (E). At A; = ±fcc, £b coincides with ±A2. 
For nonzero (3, Sb becomes exactly zero for \ky \ > ki as 
shown in curves a in Figs. 1(D) and (E). The present line 
shapes of are completely different from those of s -I- p- 
wave superconductors. On the other hand, for A^ > At, 
Eb = Q for any ky similar to the case of spin-singlet d^y 
or spin-triplet p^^-wave pairing |16l. Il7|. 

We notice here that the zero energy bound state for 
^ith fc+(2. = fc-(2. = . /fc2 _ fc2 for Ifcj^l < fci(2) I fcy I > A:i is a Majorana bound state. The wave function 
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and fcj 
and fc 



~ ^1(2) fo'^ l^yl > ^1(2)' 
(±fc^2)3;j fcy)- Here, ki and k2 are the 
Fermi wavenumbers for the smaller and larger Fermi 
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surface given by —mX/h^ + J {m\/Ti^Y + 2m/i/?i^ and 



mX/lT? + y {mX/TT?Y + 2m^/h?, respectively. The wave 
functions are given by '^V'l^ — (ij iaj~£ri±, Fij-) 

and "^"02^ — (l; *ck2±' ^o^^trai, — r2±) with 
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1(2)± 
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and Q!]^(2)± — (^^1(2)2; ^*^i/)/^i(2)- E is the quasiparticle 
energy measured from the Fermi energy. 

Postulating vE's(2;) = at x = 0, we can determine the 
ABS. We consider the case for \ky\ < k2. We first focus 
on the ABS for d^y + p-wave case. For Af > A^, the 
dispersion Eb of ABS is given by 



£b = 



(fcT+fe^TfeTi; '^'^ ^ ''^y' - (4) 

ki < \ky\ 



with 7 = (fci/fc2 + fc2/fci) + (As/At)(A:2/fci - 
/C1//C2), V = [At(l - h/k2) + A,(l + 

fcl/fc2)]/{ At[l + (fci/fc2)^] + A,[l - (fci/fc2)^]} , 

kc = ki^Atil - k,/k2) + AJ^At + A,[l ~ {ki/k2Y]. 
On the other hand, for A^ > Af , the resulting £b is given 
by Eb — 0. The dispersion Eb of ABS changes drastically 
at As = Af, where one of the energy gaps, i.e. A2, 
becomes zero. It should be remarked that the present 
ABSs do not break the time reversal symmetry. 

The resulting Eb is plotted for various cases in Fig. 
1 with Ag = As -|- Af. For convenience, we intro- 



duce dimensionless constant {3 = 2mX/{h kf) with k 



2mfi/h^. We also plot Ai and A2. Both Ai and A2 

become zero at fcj, = 0. At j/cyj = fc2, A2 is always zero. 
However, Ai then becomes zero only for (3 = 0. First, 



for the zero energy edge state "i^miky) can be written as 
'^^m{ky) = {ui{ky),U2{ky),vi{ky),V2(ky)) whcic 



ui{ky) = -i(jV2{ky) = 
U2{ky) = iavi{ky) 



{afi ~ /?i/2) exp(zfc^ y-if) 



aa 

ifi +/32/2)exp(^fcy^/- zf ) 



^k^-kD/ki, (3i = {aky/k2 + 1), h 



with a = {ky 

{a + ky/k2) and a = sgn{ky). The functions /i and 
/2 decays exponentially as a function of x and are even 
(3) function of ky. The Bogoliubov quasiparticle creation 
operator for this state is constructed in the usual way as 

J^ky) = Ul(ky)c\{ky) + M2(fcy)c|(fcy) -|- Vl {ky)c^ {- ky) + 

V2{ky)ci{—ky). Since ui{ky) = vl(—ky) and U2{ky) = 
V2{—ky) are satisfied, it is possible to verify that ^^ky) = 
l{~^ky). This means the generation of Majorana bound 
state at the edge for \ky\ > ki. For As > Af, a similar 
Majorana bound state also appears for \ky\ > ki. On 
the other hand, for \ky \ < ki, Majorana bound state has 
double branches and it is reduced to be conventional zero 
energy ABS. 

Unlike Majorana fcrmions studied before 



m m, the 



present single Majorana bound state is realized with time 
reversal symmetry. The TRI Majorana bound state has 
the following three characteristics, a) It has a unique 
flat dispersion: To be consistent with the time-reversal 
invariance, the single branch of zero mode should be sym- 
metric under ky — ^ ~ky. Therefore, by taking into ac- 
count the particle-hole symmetry as well, the fiat disper- 
sion is required. On the other hand, the conventional 
time-reversal breaking Majorana bound state has a lin- 
ear dispersion, b) The spin-orbit coupling is necessary 
to obtain the TRI Majorana bound state. Without spin- 
orbit coupling, the TRI Majorana bound state vanishes, 
c) The TRI Majorana bound state is topologically stable 
under small deformations of the Hamiltonian ([1]). 

We also calculate ABS for rfj.2_,y2 + /-wave case. In 
this case, ABS exists only for As < Af. In Fig. 2, Eb is 
plotted similarly to Fig. 1. As a reference, corresponding 
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FIG. 1: (Color online) Andreev bound state eb, effective pair 
potentials for each Fermi surface Ai and A2 are plotted for 
dxy +p-wave case as a function of fey/fe- P — for panels A, 
B and C. /3 = 0.5 for panels D, E and F. At = Aq, A^ = for 
A and D. At = O.8A0, A^ = 0.2Ao for B and E. At = 0.2Ao, 
As = O.8A0, for C and F. In all panels, curves a (solid line), 
b (dotted line) and c (dashed line) denote £i,/Ao, Ai/Ao and 
A2/A0, respectively. 



Eb is also shown for s + p-wave case. Helical edge modes 
around ky = exist and Sb is absorbed into continuum 
levels for \ky\ > ki. These features are similar to those 
of s + p-wave case. However, the number of crossing 
points of Eb is five for (1^2 _y2 + /-wave case reflecting the 
complex fc-dependence of the pair potential. The overall 
line shapes of st (curve a) in Fig. 2(A) is significantly 
different from corresponding Sb (curve a) in Fig. 2(B). 
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s+p -wave 
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FIG. 2: (Color online) Similar plots to Fig. 1 with (3 = 0.5 for 
dj;2_y2 +/-wave (A) and s-l-p-wave case (B) with At — O.8A0, 
As = O.2A0. In all panels, curves a(solid line), b(dotted 
line) and c(dashed line) denote et/Ao, Ai/Ao and A2/A0, 
respectively. 



It is very interesting to clarify how the above novel 
types of ABS are reflected in the charge transport prop- 
erty . The Hamiltonian in a normal metal is given 
by putting A(fc) = and A = in Hs. We assume an 
insulating barrier at a; = expressed by a delta-function 
potential US{x). The wave function for spin 7 = (t,i) in 



the normal metal ^^ix) is given by 
'ifN{x)^cxp{ikFyy)[{ipif+ ^ a^,pVap) exp(ifc_Fa;a;) 



^ b^,p'4jbpexp{-ikFxx)] 



(7) 



p=t,i 



with ^V»t= =(1,0,0,0), ^V^4=^V'^,i =(0,1,0,0), 

^Vat =(0,0,1,0), and ^ijjai =(0,0,0,1). The corre- 
sponding ^'s(a;) is given by Eq. ([2|). The coefficients 
a^y^p and bj^p are determined by the boundary condi- 
tion *Ar(0) = *s(0), and hvsx'^siO) - hvNx'^NiO) = 
— 2«[/f3^s(0) with fivg^jy-f^ — dHs(^iq)/dkx, and diago- 
nal matrix given by fs — diag(l, 1, —1, —1). 

The quantity of interest is the angle averaged charge 
conductance ac given by 



Ili%fNc{<P)d(p 



/c(0) = [2 + V(|a,,,p-|6^,p 



(8) 



(9) 



where fNc{4') denotes the angle resolved charge conduc- 
tance in the normal state with A(fc) = 0. Here, (j) denotes 
the injection angle measured from the normal to the in- 
terface with sm(j) = ky/kf. To characterize transparency 
of the junction interface, we introduce dimensionless con- 
stant Z = 2mU/h^kf. 

We plot bias voltage eV — E dependence of (Jc for 
dxy + p-wave case in Fig. 3 for various Z . First we 
concentrate on low transparent junction with Z — hy 
changing the value of Ag and A^. At A( = A^, one of 
the energy gap of the Fermi surface closes corresponding 
to the quantum phase transition. Then, the resulting ac 
has a gradual change from the quantum critical point. 
For the case without spin-orbit coupling (/3 = 0) with 
At > As, (JC has a gap like structure around zero bias 
due to the absence of Majorana bound state as shown 
in Figs. 1(A) and 1(B). For > At, ZBCP appears 
reflecting the zero energy ABS 17[. In the presence of 



spin-orbit coupling, uc always has a ZBCP independent 
of the ratio of A^ and At as shown in Fig. 3(B). For 
At > Ag, the ZBCP originates from purely TRI Majoana 
bound state. The width of the ZBCP for At > A^. is 
enhanced with the increase of /3, since the region of ky 
where the TRI Majorana bound state exists is expanded 
with p. For Ag > At, both the conventional ABS and 
TRI Majorana bound state contribute to the formation of 
ZBCP. We also plot corresponding ac for high {Z ^ 1) 
and intermediate {Z = 2) transparent junctions. For 
At > As , ac has a broad dip- like structure around eV — 
for Z = 1, while it is slightly enhanced around eV = 
for Z = 2 (curves a and b in Figs. 3(C) and (D)). On 
the other hand, for A^ > At, ac always has a ZBCP 
(curves d and e in Figs. 3(C) and (D)). The presence of 
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FIG. 3: (Color online) Tunneling conductance ac for dxy +p- 
wave case. A: Z = 5, ^ = 0, B: Z = 5, /3 = 0.5, C: Z = 1, 
P = 0.5, D: Z = 2, /? = 0.5. a(solid line): A* = Aq, A, = 0, 
b(dotted line): At = 0.8 Ao, A^ = 0.2Ao, c(dashed line): 
At = 0.5Ao, As = 0.5Ao, d(tliin solid line): At = 0.4Ao, 
As = O.6A0, and e(dot-dashed line): At = 0.2Ao, A^ = 
O.8A0. 
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FIG. 4: (Color online) Tunneling conductance gc with /3 = 
0.5 and Z = 5. A: d^2_y2 +/-wave and B: s+p-wave. a(solid 
line): At = Ao, A^ = 0, b(dotted line): At = O.8A0, A^ = 
O.2A0, c(daslied line): At = 0.5Ao, As = 0.5Ao, d(tliin solid 
line): At = 0.4Ao, As = O.6A0, and e(dot-dashed line): At = 
O.2A0, As = O.8A0. 



TRI Majorana bound state gives a clear ZBCP with the 
increase of Z. As a reference, the tunnehng conductance 
ac for d^2_y2 +/-wave and s+p-wave cases are plotted in 
Fig. 4 for Z = 5. ABS exists only for < Aj. The ac 
for dx2-y2 + /-wave has a ZBCP splitting reflecting the 
complex dispersion Sb shown in Fig. 4(A). On the other 
hand, for s + p-wave case, ac has a broad ZBCP shown 
in Fig. 4(B). Summarizing Figs. 3 and 4, ac for each 
paring state are qualitatively different from each other, 
which can be used to identify these pairings. 

In conclusion, we have studied the ABS and resulting 



charge transport for dxy +p-wave and dx2-y2 + /-wave 
superconductors. We find that the obtained dispersion 
of ABS in both cases have an anomalous structure. For 
dxy +p-wave case, a novel TRI Majorana bound state is 
generated due to the spin-orbit coupling. The resulting 
charge conductance can serve as a guide to identify the 
TRI Majorana bound state and paring symmetry of NCS 
superconductors by tunneling spectroscopy. 
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